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Recent results on two-loop calculations in the MSSM are reviewed. The computa- 
tion of the QCD corrections to Ap, Ar and to the mass of the lightest Higgs boson 
in the MSSM are presented. 

1 Introduction 

1.1 Motivation 

Supersymmetry ISUSY) , in particular the Minimal Supersymmetric Standard 
Model (MSSM), El is widely considered as the theoretically most appealing ex- 
tension of the Standard Model (SM). It predicts the existence of scalar partners 
Jl,Jr to each SM chiral fermion, and spin-1/2 partners to the gauge bosons 
and to the scalar Higgs bosons. So far, the direct search for SUSY particles 
has not been, successful. One can only set lower bounds of 0(100) GeV on 
their masses El. 

An alternative way to probe the MSSM is to search for virtual effects of the 
additional particles by means of precisian observables. Here most of the one- 
loop calculations are already available Era. In order to reach the same standard 
of theoretical predictions as in the SM, two-loop calculations are required. 

Another way to test the MSSM is to search for the lightest Higgs boson, 
which is at tree-level predicted to. be lighter than the Z boson. The one-loop 
corrections are known to be hugeEl. Hence a two-loop calculation is inevitable 
to have the mass predictions under control. 

Since the dominant two-loop contributions for electroweak precision ob- 
servables as well as for the Higgs sector are expected to originate from the 
QCD corrections, we concentrate on 0(a s ) calculations in this article. 

1.2 Notations 

Contrary to the SM,_in the MSSM two Higgs doublets are needed. The Higgs 



potential is given byl 




(1) 




(HxHi — H2H2) 2 + — \HiH2\ 2 , 
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where m^, TO12 are soft SUSY-breaking terms, and g, g' are the SU (2) and U(l) 
gauge couplings. 

The doublet fields Hi and H 2 are decomposed in the following way: 
R _ (H{\ _ (v 1 + ($f + ixf)/V2 

The vacuum expectation values V\ and v 2 define the angle f3 via 

tan^ = — . (3) 

Vl 

In order to obtain the CP-even neutral mass eigenstates, the rotation 



H°\ _ f cos a sin a \ / 



h° ) \ — sin a cos a / \ $2 



(4) 



is performed. The potential (gj) can be described with the help of two inde- 
pendent parameters: tan/? and M\ — — m 2 2 (tan/3 +cot/3 ), where Ma is the 
mass of the CP-odd A boson. 

SUSY associates a left- and a right-handed scalar partner to each SM 
quark. The current eigenstates, ql and qu, mix to give the mass eigenstates 
qi and q 2 . In the MSSM the q mass eigenstates and their mixing angles are 
determined by diagonalizing the following mass matrix (e q and If are the 
electric charge and the weak isospin of the partner quark, Mq L , Mq R and A q 
are soft SUSY-breaking terms) 

2 / Ml + m 2 + cos 2/3(If - e,^) Mf m g M^ R \ 

(5) 

with s 2 ^ = 1 — c 2 ^ = sin 2 9w ■ Especially for the i — b sector the parameters 
are given by M t Li? = A t — /i cot/3 and M^ R = Ab — /itan/3 . Furthermore, 
SU(2) gauge invariance requires M^ L — My L at the tree-level. 

Due to the large value of the top-quark mass mt, the mixing between the 
left- and right-handed top squarks fx and tn can be very large. 



2 Evaluation of the two loop diagrams 

For precision observables like the p parameter and Ar, which fixes the higher 
order relation between My/, Mz, Gf and a, the leading contributions can be 
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expressed in terms of self-energies. The corrections to the Higgs masses are 
given partly in terms of self-energies and partly in terms of tadpole diagrams. 
Hence for the corrections presented in this paper, it is sufficient to concentrate 
on one- and two-point functions. 

One problem in dealing with two-loop calculations is the large number of 
topologies and the even much larger number of diagrams. Thus we make exten- 
sive use of computer algebra programs: The two-loop diagrams for a specific 
process, including also the required jcounterterm diagrams, are generated with 
the Mathematica package FeynArtsu. The model file we programmed for this 
purpose contains, besides the SM propagators and vertices, the relevant part 
of the MSSM Lagrangian, i.e. all SUSY propagators (t i , t 2 , b\ , b 2 , g ) needed for 
the QCD-corrections and the appropriate vertices (gauge boson-squark ver- 
tices, squark-gluon and squark-gluino vertices). FeynArts inserts propagators 
and vertices into the graphs in all possible ways and creates the amplitudes 
including all symmetry factors. 

As shown inB, an arbitrary two-loop self-energy diagram can be reduced 
to a set of scalar one- and two-loop integrals by means of two-lojQp tensor 
integral decomposition. Besides the well known functions A« and BqB : this set 
consists of the functions T12345, T1234, T234, T134, defined asQ 



Ti2...5(k 2 , mi,m 2 , . . . , m 5 ) = J d D q x d D q 2 -. 
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■[kl- m \][kl- m l]...[ki-rniY 

(6) 

with hi as internal momenta, as internal masses and q\ 2 as integration 
momenta. The decomposition of the two-loop diagrams and counterterms is 
performed with the Mathematica package TwoCalcu. 

The scalar integrals can be split into a divergent and a finite part, perform- 
ing an expansion in 5 — (4 — D)/2 (where D is the dimension of space-time). 
This allows, after including the proper counterterm diagrams, an algebraical 
check of the finiteness. The finite result can now be transformed into a FOR- 
TRAN code, which makes a fast numerical evaluation possible. 



3 Dimensional Reduction (DRED) 

3.1 What is the problem? 

In SUSY theories the number of bosonic degrees of freedom is equal to the 
number of the fermionic ones. If the regularization is performed through di- 
mensional regularization (DREG), all the fields are treated D-dimensional. A 
problem arises through the fact that a vector in D dimensions has a different 
number of d.o.f. than a vector in 4 dimensions, thus DREG breaks SUSY. 
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As a solution DRED was proposed EJ. Here the indices of the fields (and 
corresponding matrices) are kept 4-dimensional, the integrals and momenta 
are treated D-dimensional. Alternatively, it is possible to use the so called 
e-scalarsO. 

For the program TwoCalc we implemented the option to perform the cal- 
culation either in DREG or in DRED. 

3.2 What is changing? 

As explained in the previous chapter, a diagram is decomposed into a set of 
basic scalar integrals. For the coefficients of these integrals the calculation in 
DRED can give an extra contribution ~ 5 with respect to DREG. The changes 
for a diagram are the following: 

• ^—divergencies: no change from DREG to DRED. 

• ^-divergencies: for one-loop diagrams no change occurs, for two-loop 
diagrams, however, there is an extra contribution, i.e. it is possible that 
a calculation is divergent in DREG but finite in DRED (see Sec. [?]). 

• finite part: for one loop diagrams as well as for two-loop diagrams there 
arises an extra contribution. Thus it is possible that a one-loop calcula- 
tion is finite in DREG and in DRED, but has different results. 

4 Renormalization 

For our two-loop calculations for precision observables, renormalization is re- 
quired up to the two-loop level. 

• Gauge boson sector: Two-loop renormalization is required. For the W- 
boson we have 

mass renormalization: — > Myy + 8Mw\ + SM^ 2 

field renormalization: — > z\Jj 2 W^ (7) 

with z\l 2 = (1 + SZw.i + SZ W ,2) 1/2 = 1 + T.SZ w ,i - + t;& z w,2 

A o ' Z 

In the on-shell renormalization scheme one has: 

*Af^ il = ReSw,i(M^), 

SM W2 = Re£vK,2(-Mw) + combination of one-loop terms, 
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6Z Wtl = -ReZ' Wtl (M^), (£'(fc 2 ) := — E(i 2 )) (8) 
&Zw,2 = — R-eS' H / 2 (-^H / ) + combination of one-loop terms. 
For the Z this is performed analogously. 

• Quark sector: One-loop renormalization as in the SM has to be inserted. 

• Squark sector: One-loop renormalization is sufficient. It can be per- 
formed in terms of masses and mixing angle: 

TO? -> TO? + <5to? , e~ q 9g + 56 q. (9) 

For the masses on-shell conditions are used. The renormalized mixing 
angle can be defined by requiring that the renormalized squark mixing 
self-energy T, T ^~ 2 (q 2 ) vanishes when the q% is on-shell: 

£m| = Ej.(m|), 69g= 2 - ^ 2 Sft^mfJ. (10) 

91 92 

However, one has to keep in mind that there are relations between the 
soft SUSY-breaking terms in the squark mass matxices, which also have 
to be renormalized for a two-loop calculation (see til). 



5 QCD corrections to Ap 

5.1 Definition and one-loop results 

The deviation of the p parameter from unity parametrizes the leading universal 
corrections induced by heavy fields in electroweak amplitudes, caused by a 
mass splitting between the partners in an isospin doublet. Compared to this 
correction all additional contributions are suppressed. Ap gives an important 
contribution to electroweak observables via 



A(Ar)«-4 1 Ap, AM W ^^f ^L_ Ap, (111 



Asin 2 flff«- f w *\ Ap. (12) 
c w s w 

In terms of the transverse parts of the W- and Z-boson self-energies at 
zero momentum-transfer, the p parameter is given by 

i a Mo) *V(o) 

p = l^' Ap = ^if-^qr- (13) 
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Inserting the MSSM contribution to the W- and Z-boson self-energies 
(diagrams of Fig. [l]) and neglecting the mixing in the sbottom sector, one 
obtains for the contribution of the t/b doublet at one-loop order (only 6l 
contributes for 0^ — 0) : 



sin 2 Of cos 2 Of F (ml , mf ) 



+ cos 2 t -Fo(m 2 i ,m? i ) + sin 2 ^ F (m 2 2 , m? j] , (14) 

with Fq(x, y) = x + y In—, (15) 

x-y y 

where Fo has the properties Fq(x,x) — 0,Fo(a;,0) = x. Therefore, the con- 
tribution of a squark doublet becomes in principle very large when the mass 
splitting between the squarks is large. This is exactly the same situation as in 
the case of the SM where the top/bottom contribution to the p parameter at 
one-loop order readsll3 

Ap ° M = ^fe F ° (m? ' W ' ) - (16) 

For m t 3> m this leads to the well-known quadratic correction Apg = 
3G F m 2 t /(8V2ir 2 ). 



iAA/WW) WVWW 

V '\ / V 

aaaa/wwMaAaaavvwv- 



Figure 1: Feynman diagrams for the contribution of scalar quark loops to the gauge boson 
self— energies at one-loop order. 

In Fig. U we display the one-loop correction to the p parameter induced by 
the t/b isodoublet. Due to the SU(2) symmetry and for the sake of simplicity 
we concentrated on the case Mf — Mf — Mr — Mr = mg. In this scenario, 
the scalar top mixing angle is either very small, 9f ~ 0, or almost maximal, 
Of ~ — 7r/4, in most of the MSSM parameter space. The contribution ApQ USY 
is shown as a function of the common squark mass rriq for tan/3 = 1.6 for 
the two cases M^ R = (no mixing) and M^ R — 200 GeV (maximal mixing); 
the bottom mass and therefore the mixing in the sbottom sector are neglected 
leading to = ~ rriq. Here and in all the numerical calculations in 
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this paper we use m t = 175 GeV, M z = 91.187 GeV, M w = 80.33 GeV, and 
a s = 0.12. The electroweak mixing angle is defined from the ratio of the vector 
boson masses: Syr = 1 — M^/M§. 




5.2 Two-loop results 

The QCD corrections to Ap via the squark contributions to the vector boson 
self-energies, Fig. |j| have been derived recently E3'L3. The diagrams can be 
divided into three different classes: the pure scalar diagrams (Fig. |^a, 2. line), 
the gluon-exchange diagrams (Fig. |3|a, 1. line), and the gluino-exchange dia- 
grams (Fig. |^b). These diagrams have to be supplemented by counterterms for 
the squark and quark mass renormalization as well as for the renormalization 
of the squark mixing angle. 

The three different sets of contributions together with their respective 
counterterms are separately gauge-invariant and ultraviolet finite. For the 
gluon-exchange contribution we have only considered the squark loops, since 
the gluon-exchange in quark loops is_jiist the SM contribution, yielding the 
result Apf M = -Apg M |^(l + 7r 2 /3)H 

The pure scalar diagrams give zero contribution: partly they only con- 
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Figure 3: Typical Feynman diagrams for the contribution of scalar quarks and gluinos to 
the W/Z-boson self— energies at the two— loop level. 



tribute to the longitudinal part of the vector boson self-energies, or they cancel 
exactly with their corresponding counterterm diagram. 

The gluon-exchange contribution can be cast into a very simple formula: 



ApfUSY 



G 



F a s 



gluon 



4\/27r 3 

i „„„2 



sin Or COS 



Tflr , mr 
12 ' b L 



(17) 



with Fi{x,y) = x + y — 2- 



In* 
2/ 



2 + - In - 

V V. 



(x + y)x 2 2 x 



(x-y) 2 



In - - 2(x 



j/)Li 2 (l--) 

y 



(18) 



where F\ has the properties F\(x, x) = 0, F\(x, 0) = x(l + 7r 2 /3). The two-loop 
gluonic SUSY contribution to Ap is shown in Fig. |] for the same parameters 
an in Fig. 0. 
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Like for the gluon-exchange contribution we have also derived a complete 
analytic result for the gluino-exchange contribution. This is, however, very 
lengthy and we therefore present here a numerical result in Fig. |5| for the same 
parameters as in Fig. [|, both for a light and a heavy gluino. (The complete 
expression is available in Fortran and Mathematica format from the author.) 

The complete result has been derived in DREG and DRED. Both methods 
yield the same result, since the only diagram which differs in both schemes is 
the gluon-exchange in the top self-energy. This, however, contributes to the 
SM part only. 

To summarize, the gluonic contribution amounts to ~ 10 — 15% of the 
one-loop result. Contrary to the SM two-loop corrections, they have the same 
sign, resulting in an enhancement in the sensitivity to virtual squark effects. 
The gluino-exchange contribution is in general smaller compared to the gluon 
result and has different signs. It competes with the gluon contribution only for 
rrig and close to their lower experimental bound. There both correction 
add up to ~ 35% for maximal mixing. The gluino contribution decouples for 
heavy m g . We confirmed this by an explicit series expansion in rrig starting 
with an coefficient of 0(l/m g ). 
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Figure 5: Contribution of the gluino exchange diagrams to Ap^USY^ f or t wo values of nig 
in the scenarios of Fig. M. 



6 QCD corrections to Ar 

Ar fixes the higher order relation between Mw, Mz, Gf and a: 



M w 



1 



Ml 



ira 
V2G F 



(1 + Ar), 



(19) 



where all the radiative corrections are contained in Ar. At two-loop in leading 
order the SUSY contribution is given by 



Ar SUSY 



= £i(0)-T 



5M% 



E W (0)-5M^ | 2 c w S 7Z (0) 
s w Mf : 
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(20) 

with SM^ = ReSy(M^) if on-shell renormalization is used. The SUSY QCD 
corrections thus only enter via the self-energies. 

Since Ar determines the prediction for the mass of the W-boson, it is 
crucial to know Ar with high precision. The complete MSSM one-loop result 
is available 0. But at least through Ap (|l2|) also a non-negligible two-loop 
correction is expected. 

The diagrams involved here are the same as in the Ap case but in gen- 
eral have non vanishing external momentum. Therefore they have a much 
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more complicated structure. As a first step an analytical result for the gluon- 
cxchange contribution has been derived. Due to the fact that nig ^ 0, the 
gluino contribution can not be expressed analytically. Since for large nig the 
gluino decouples, the gluon-exchange contribution is a good approximation for 
the case of a heavy gluino. 

The result for a general seLtenergy Ey(fc 2 ) at the two-loop level can be 
cast into a very compact formli3, consisting of the scalar one-loop functions 
Aq,Bq as well as of the two-loonJunctions T11234', J1234', Ti23'45 and T234'- 
The T-functions can be found inllH. After expanding the complete result for 
Ar SUSY in 6, we found a finite result for the gluon-exchange contribution. 

The result is presented in Fig. ^ for the same parameters as in Fig. |[ 
including also the Ap-approximation defined in (|l2|). 



CO 
CO 




— 0.0020 



100 



1.6 

OGeV 

GeV (Ap-approx.) 
200 GeV 

200 GeV (Ap-approx.) 



200 



300 

m~ 



400 



500 



Figure 6: Contribution of the gluon-exchange diagrams to Ar s 
approximation given for the scenarios of Fig. ^. 



and the Ap- 



The two-loop correction amounts ~ 10 — 15% of the one-loop result. The 
Ap-approximation reproduces the two-loop contribution up to 10— 15%. From 
this it can be conjectured that also the gluino contribution can be well approx- 
imated with the help of Ap. Expressed in terms of the W-mass, the two-loop 
MSSM contribution corresponds to a shift in Myy of about 10 — 20 MeV for 
light squarks. 
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7 Mass of the lightest MSSM Higgs boson: QCD corrections 



7.1 Renormalization 

The MSSM predicts that one of the Higgs bosons is relatively light. At tree- 
level mh < Mz holds. However, this is changed by radiative corrections. After 
taking the one-loop correction into account, nth > Mz is possible, but a new 
upper bound is set, rrih < 15Q_GeV. At the two-loop level the upper bound is 
decreascd_t£> about 130 GeVEj. Also an RGE calculation exists at the two- 
loop level E3 which shows that is considerably reduced with respect to the 
one-loop result. 

A diagrammatic calculation has been performed to determine mh at the 
two-loop level for general values of tan/3 and Ma- In order to simplify our 
calculation, we work with the unrotated fields $i and $2 (but we obtain an 
a e g at two-loop in the end) . We only consider the t — t-sector with no mixing, 
so we set M t = 0. The leading terms are expected to originate for the <1> 
self-energies evaluated for zero external momentum (which is then similar to 
the effective potential approach). We only consider the Yukawa part of the 
theory, since due to the high top mass the Yukawa couplings are relatively 
large. Thus we set the gauge couplings to zero. 

The counterterm for a <& self-energy arises only from the Higgs poten- 
tial (Q) where we used on-shell renormalization for the A-boson self-energy. 
The tadpole renormalization has been chosen to cancel the tadpole contribu- 
tions, this leads to 

E$(0) = S$(0) — <5V$, (S : renormalized self-energy) 

SV* = h(Sti,6t 2 ,SMl) (21) 
with SU = -T h 5M 2 A = S A (0), 

where T,; denotes the tadpole contribution, Sti is the corresponding counter 
term. 

The QCD corrections to via the squark contributions to the $i self- 
enecgies, Fig. |?] (we omitted the tadpole diagrams), are worked out in detail 
in El. Again the diagrams can be divided into three different classes: the 
pure scalar diagrams (Fig. |^a), the gluon exchange diagrams (Fig. |^b), where 
also the quark loops have to be taken into account, and finally the gluino- 
cxchange diagrams (Fig. |c). These diagrams have to be supplemented by 
counterterm insertions for the squark and quark mass renormalization as well 
as for a renormalization of the squark mixing angle. 

Contrary to the Ap case the three sets of diagrams do not independently 
give an ultraviolet finite result for the renormalized $ self-energy. In addition it 
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Figure 7: Typical Feynman diagrams for the contribution of quarks and scalar quarks with 
gluons and gluinos to Higgs-boson self-energies at the two-loop level. 



is necessary to make use of DRED. By explicit calculations both in DREG and 
in DRED, we found a finite result only for the DRED. The same observation 
has been made also inlia for the two-loop calculation of the upper bound for 
m h . 

7.2 Calculation of 

The mass matrix at two-loop order is given in terms of the tree-level masses 
and of the one- and two-loop renormalized $ self-energies which can be diag- 
onalized by the angle a e g 



Hi = 



- y; 1 - 

^$1*2 



- v 2 



'*1*2 



yi _ y2 



)»7i 



Y 1 



Y 2 

^$2 



'H 







where m^.h denote the heavy and light Higgs masses at two-loop. 

Due to the gluino-exchange contribution the analytical result is very com- 
plicated, it will be given inE3. The numerical results are displayed in Fig. || for 
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raj, at the tree-, the one- and the two-loop level as a function of rriq for the 
same scenario as in Sec. [3] but with tan/3 = 1.6 and 40, rrig = 500 GeV and 
M A = 200 GeV (M t LR = GeV.) 



120 - 



100 



tan [3 = 40 



> 80 : 

CD 
O 

— 60 - 



40 



20 



tanp = 1.6 



m 5 = 500 GeV, M A = 200 GeV 

tree level 

1 loop 

— 2-ioop ; 



100 200 300 400 500 600 700 800 900 1000 

m~ [GeV] 



Figure 8: The mass of the lightest MSSM Higgs boson rrth displayed as a function of rriq 
for two values of tan j3 , m g = 500 GeV , M A = 200 GeV 



The two-loop result lowers the one-loop mass by about 15 GeV. The 
decrease is even somewhat larger for small tan (3 . This strengthens for the 
small tan j3 scenario, which is favoured by SU(5) GUT scenarios, the discovery 
potential of LEP2. 
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